Half-metallic ferromagnetism in binary compounds 
of alkali metals with nitrogen: Ab initio calculations 
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The first-principles full-potential linearized augmented plane-wave method based on density 
functional theory is used to investigate electronic structure and magnetic properties of hypothetical 
binary compounds of l A subgroup elements with nitrogen (LiN, NaN, KN and RbN) in assumed 
three types of cristalline structure (rock salt, wurtzite and zinc-blende). We find that, due to the 
spin polarized p orbitals of N, all four compounds are half-metallic ferromagnets with wide energy 
bandgaps (up to 2.0 eV). The calculated total magnetic moment in all investigated compounds 
for all three types of crystal structure is exactly 2.00 /j,b per formula unit. The predicted 
half-metallicity is robust with respect to lattice-constant contraction. In all the cases ferromagnetic 
phase is energetically favored with respect to the paramagnetic one. The mechanism leading to 
half-metallic ferromagnetism and synthesis possibilities are discussed. 



PACS numbers: 71.15.Mb, 71.20.Dg, 72.25.Ba, 75.50.Cc 
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I. INTRODUCTION 

Half-metallic (HM) ferromagnets are materials in 
which, due to the ferromagnetic decoupling, one of the 
spin subbands is metallic, whereas the Fermi level falls 
into a gap of the other subband. The concept of HM 
ferromagnet was first introduced by de Groot et al. in 
1983 [l| on the basis of band structure calculations for 
NiMnSb and PtMnSb semi-Heusler alloys. HM ferromag- 
nets are considered as promising materials to exploit the 
spin of charge carriers in new generations of transistors 
and other integrated spintronic devices @, in particu- 
lar, as a source of spin-polarized carriers injected into 
semiconductors, since only one spin channel is active dur- 
ing charge transport, thus leading to 100% spin-polarized 
electric current. 

Since 1983 many HM ferromagnets have been theoret- 
ically predicted, but very few of them found experimen- 
tal confirmation like metallic oxides Cr02 [3J] and Fe3C>4 
[U or manganese perovskite La .7Sr .3MnO3 @. Many 
HM ferromagnetic materials were predicted in transition- 
metals pnictides [a, 0] and chalcogenides [g, @| by means 
of first-principles calculations. 

Recently, an unusual class of ferromagnetic materials 
(lCi UlL [l2l . [la . LLJ, LL5[ , which do not contain transition- 
metal or rare-earth atoms, has been predicted and anal- 
ysed theoretically by ab initio calculations. In Ref. [l(| 
the authors present ab initio calculations for CaAs in 
zinc-blende structure, where magnetic order is created 
with the main contribution of the anion p electrons ( "p- 
electron" ferromagnetism). More comprehensive study 
was made in , where the authors investigate p-electron 
ferromagnetism in a number of tetrahedrally coordinated 
binary compounds of I/II-V elements. The characteris- 
tic feature of this class of materials is the integer value of 
the total magnetic moment per formula unit, which, in 
some combinations of elements, can be as large as 3 fiB 
15]. Results presented in Ref. [ll|, where magnetic and 



structural properties of II^-V nitrides have been inves- 
tigated using ab initio methods, motivated us to check 
possibility of finding half-metallic ferromagnetism in \ A - 
N binary compounds. Another motivation was that the 
atomic bonds in I^-N nitrides are supposed to be mostly 
ionic in nature (due to significant difference in electroneg- 
ativity between I atoms and nitrogen atom) which was 
essential for appearing of p-electron ferromagnetism in 
all previously considered HM cases 

In this paper we present electronic structure and mag- 
netic properties of hypothetical 1 A -N binary compounds 
(LiN, Nan, KN and RbN) with the rock salt (RS), 
wurtzite (WZ) and zinc-blende (ZB) crystalline struc- 
ture, calculated by means of first-principles full-potential 
linearized augmented plane-wave method. We find that 
all four compounds in all three types of structures are HM 
ferromagnets with robust half-metallicity against lattice 
compresion in the range from 2% for LiN (RS) up to 50% 
for NaN (WZ). This is crucial parameter for the practice 
of epitaxial growth, e.g. by means of MBE or MOCVD, 
to fabricate magnetic ultrathin layer structures for spin 
injection into suitable semiconductor substrate. 

The paper is organized as follows. Section II shows de- 
tails of our calculation method, while section III presents 
calculated total energy and band structure. In section 
IV we analyse density of states and the origin of half- 
metallic ferromagnetism. The paper is concluded with 
section V. 



II. COMPUTATIONAL METHOD 

All calculations were performed using WIEN2k code 
16] which implements the full-potential linearized aug- 
mented plane wave (FLAPW) method [l7j]. Exchange 
and correlation were treated in the local spin density ap- 
proximation (LSDA) by adding gradient terms. This 
GGA approximation was used in the Perdew-Burke- 
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Ernzerhoff [18] parametrization. It should be mentioned 
that LSDA gives the magnetic moments in a very good 
agreement with experiment but underestimates the lat- 
tice constants in the case of transition metals. On the 
other hand, gradient corrections significantly reduce this 
error and give the correct phase stability but tend to 
overestimate the magnetic moment. We checked this de- 
ficiency of GGA in our calculations and it turned out 
that it has a small efect on the considered systems. 

The convergence of the basis set was controlled by the 
cut-off parameter RK maK = 8 together with 5000 k-point 
mesh for the integration over the Brillouin zone. The an- 
gular momentum expansion up to I = 10 and G max =12 
a.u. -1 (G max =14 a.u. -1 in case of RbN) for the poten- 
tial and charge density was employed in the calculations. 
Self-consistency was considered to be achieved when the 
total energy difference between succeeding iterations is 
less then 10~ 5 Ry per formula unit. Geometry optimiza- 
tion was performed allowing all atoms in the unit cell to 
relax, constrained to the initially assumed crystal sym- 
metry. 



III. TOTAL ENERGY AND BAND STRUCTURE 

For all four compounds investigated by us the total en- 
ergies E have been calculated as a function of the lattice 
parameter a for three crystal structures, namely RS, WZ 
and ZB, in ferromagnetic and paramagnetic state. We 
also optimized lattice constants a by minimization of the 
total energy E. Figure 1 presents total energy E versus 
lattice constant a for all three crystal structures. As one 
can see in all compounds the RS strucure is the most 
energetically favored, thus, indicating the ionic bonding 
as prevailed coupling mechanism. Also, in all the cases, 
the paramagnetic phase is higher in energy than the cor- 
responding ferromagnetic one. Calculated equilibrium 
lattice constants ao can be found in Table I as well as 
differences between total equilibrium energies of ferro- 
magnetic and paramagnetic states AE t '^. 

The calculated total magnetic moment in all investi- 
gated compounds for all three types of crystal structures 
is exactly 2.00 /ib per formula unit. An integer value of 
magnetic moment is a characteristic feature of HM ferro- 
magnetism [Tlj . Table I shows total values of magnetic 
moments for all compounds as well as contributions from 
l A atoms (Li, Na, K and Rb), N atom and interstitial 
region to total magnetic moment. The main contribu- 
tion in all cases comes from nitrogen anion ranging from 
1.51 ^ B for LiN in ZB to 1.85 ^ B for KN in RS and ZB, 
which confirms the general feature that in HM ferromag- 
netic compounds the most of resulting global magnetic 
moment is carried by anion electrons. To verify that 
the integer value of total magnetic moment stabilizes the 
crystal structure, we have calculated the total energy E 
as a function of the total magnetic moment ^ to t ■ Indeed, 
Figure 2 shows that the creation of magnetic moment 
equal to 2.00 /is leads to the minimization of the total 



energy in all the cases. 

/From the application point of view it is important 
to study robustness of the half-metallicity with respect 
to lattice constant. Figure 3 shows magnetic moment 
as a function of lattice constant for all four compounds 
in three crystal structures. In the case of LiN the total 
magnetic moment remains integer until the lattice con- 
stant is compressed to critical value of 2.75 A, 3.75 A and 
4.25 A for WZ, ZB and RS structures, respectively. In 
case of NaN, KN and RbN the values can be read from 
Figure 3. In all studied compounds the most promising is 
wurtzite structure, where half-metallicity is maintained 
up to contraction of the lattice parameter of 9% (LiN), 
54% (NaN), 22% (KN) and 21% (RbN). Wide bandgap 
semiconductors A1N, GaN and ZnO in WZ crystal struc- 
ture are considered as strong potential materials for spin- 
tronic applications. The spin scattering relaxation time 
of charge carriers during the transport process in these 
materials is strongly increased as compared to GaAs [li| . 
This means longer lifetime of electrons in particular spin 
state and longer spin "memory" of electrons when in- 
jected into these materials. Lattice constants of these 
wide bandgap semiconductors are close to 3.2 A, what 
is in the range of LiN and NaN HM ferromagnetism in 
WZ structure with 3% and 18% lattice mismatch, re- 
spectively. Therefore, having in mind state-of-the-art in 
producing artificial structures, at least some of the HM 
ferromagnets studied here seem to be potential candi- 
dates for epitaxial growth on A1N, GaN and ZnO sub- 
strates as layers for injection of 100 % spin polarized 
electrons. Most of ZB III-V semiconductors as well as 
narrow bandgap IV- VI semiconductors in RS structure 
(binary compounds of Pb with S, Se and Te) also sat- 
isfy lattice matching conditions as possible substrates for 
bulk growth of proposed here T*-N compounds in HM 
ferromagnetic phase. 

Figures 4 and 5 illustrate calculated spin resolved elec- 
tron band structure for RbN at equilibrium lattice con- 
stant for three types of crystal structures. In all cases 
the minority spin subband is metallic, while the major- 
ity spin subbands are separated from the Fermi level by 
wide bandgaps of 4.22 eV, 3.48 eV and 3.45 eV in RS, WZ 
and ZB crystal structure, respectively. This is in contrast 
with the 3d pnictides in which the bandgap appears in 
the minority spin subband as a consequence of d bands 
splitting [jj. Three remaining compounds (LiN, NaN, 
and KN) also demonstrate half-metallic nature with wide 
bandgaps. The HM bandgaps, separating allowed ener- 
gies for majority spin states from the Fermi level, are 
wide enough and ranging from about 0.2 eV for LiN (RS) 
to no less than 2.0 eV in the case of KN and RbN for 
three considered crystal structures. The last value is, to 
our knowledge, the largest one ever obtained for binary 
compound HM ferromagnet, thus making KN and RbN 
very promising materials for spintronic applications. The 
bands for RbN are flat and nearly dispersionless close to 
the Fermi level. The mechanism leading to this kind of 
energy bands has been discussed by other authors [ll|, [l2j 
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as an important condition for stability of HM ferromag- 
netism. 

In the case of NaN and KN the bands look very simi- 
lar to those presented in Figures 4 and 5 for RbN. For 
LiN the situation is slightly different, because the bands 
near the Fermi level are more dispersive, but far from 
breaking down HM ferromagnetism. Table II presents 
calculated values of majority spin-up E^ p and minority 
spin-down E^ n subbands main bandgaps, as well as HM 
bandgaps E^ M in the majority-spin subband. Decrease 
of main bandgaps (of direct character only in the case of 
RS structure) and increase of HM bandgaps with increas- 
ing atomic number of alkali metal elements is observed. 



IV. DENSITY OF STATES AND HM 
FERROMAGNETISM 

It is instructive to compare the band structure (Figs. 
4 and 5) with spin resolved total and partial site and 
symmetry projected density of states (DOS) presented 
in Figure 6. There are six valence electrons in RbN (Rb: 
5s 1 ; N: 2s 2 2p 3 ). Two of them occupy the low-energy 
N s states, about 11 eV below the Fermi level, and the 
remaining four electrons are mainly involved in filling an- 
ion N p states. The crystal field of cubic or hexagonal 
symmetry caused by surrounding N anions splits the Rb 
cation Ad states with three t2 g states (d xy , d yz , d zx ) lying 
lover in energy and the two e g states (d x 2 _ y 2 , d 3z 2 _ r 2) 
lying higher. Hybridization between N p states and Rb 
t2 g or e g Ad states (depending on the type of crystal 
lattice) creates both bonding and antibonding hybrid or- 
bitals. The bonding orbital, lower in energy, appears at 
the edge of normally occupied region of N p states. The 
antibonding hybrid orbital remains in the Rb ti g or e g 
manifold but pushed up in energy relative to the non- 
bonding state. 

The bands around Fermi level are formed mainly of N 
p states with small admixture of Rb d states with pre- 
dominant e g symmetry in RS structure and ti g symmetry 
for ZB and WZ structures for both spin orientations. In 
contrast to the p states, the s states are placed well below 
or high above Fermi level giving no net spin polarization. 
Similarly, the bands located at higher energies (from 4.5 
to 12.5 eV above the Fermi level) are mainly composed of 
d states belonging to t2 g and e g representations. Hence, 
the s and d states are not directly involved in creation of 
spin polarisation and ferromagnetism in RbN and other 
I^-N binary compounds. 

Existence of narrow bands near Fermi level translates 
into peaks in the DOS diagrams. According to qualita- 
tive Stoner criterion the high density of electron states at 
the Fermi level should stabilize the ferromagnetic order. 
Pressure induced lattice contraction produces changes in 
the band structure, increase of band widths with increas- 
ing kinetic energy, decrease of the spin splitting and de- 
stroying ferromagnetism. 

In the presented here results of the ground state total 



energy calculations we can see the influence of Coulomb 
intra-atomic electron repulsion. It produces energy shift 
A = 2 eV between narrow majority and minority spin en- 
ergy band states, visible as peaks at - and below Fermi 
level in the DOS diagrams, shown in Figure 6. This type 
of band structure with well separated up and down spin 
states can be attributed to the Hubbard-like interaction 
U = A/2. The parameter U = 1 eV is a measure of 
energy formation for magnetic moment localized at ni- 
trogen ion. According to Hund's rule and Pauli principle 
the limiting value of magnetic moment for p electrons is 
the value of 3 /^b- The four valence p electrons, present 
in the system, can produce magnetic moment not greater 
than (3 — 1) //b = 2 //b, what agrees with our calcula- 
tions. 

Estimation of thermodynamic stability of HM ferro- 
magnetism requires additional calculations to determine 
the strength of inter-atomic exchange coupling [20] . This 
task requires supercell calculations for planar antiferro- 
magnetic structure of the first (AF1) and second kind 
(AF2) as particular cases of spiral structure with the q 
vector [001] and 5 [111], respectively [ll|. We have per- 
formed relevant calculations only for LiN in WZ crystal 
structure. The energy differences are -Eafmi — E-pm = 
0.247 eV and £ PA ra - E FM = 0.701 eV. The estimated 
lover limit for paramagnetic Curie temperature is about 
464 K, what makes the material promising for experi- 
mental investigation. 



V. CONCLUSIONS 

On the basis of ab initio calculations employing density 
functional theory we have investigated half-metallic fer- 
romagnetism in rock salt, wurtzite and zinc-blende com- 
pounds composed of group \ A alkali metals as cations and 
nitrogen as anion. We find that, due to the spin polar- 
ized p orbitals of N, all four compounds are half-metallic 
ferromagnets with wide energy bandgaps (up to 2.0 eV). 
The calculated total magnetic moment in all investigated 
compounds is exactly 2.00 /ib per formula unit. Our cal- 
culations show that the predicted half-metallicity is ro- 
bust with respect to lattice constant contraction. The 
formation of ferromagnetic order requires large lattice 
constants, high ionicity, empty d orbitals and slight hy- 
bridization between N anion p states and I A cation d 
states with energies in vicinity of the Fermi level. It is 
interesting to note that palladium (with its eight Ad elec- 
trons) , as dopant replacing Ga atom in WZ GaN semicon- 
ductor, interacts with surrounding N atoms in a similar 
way like alkali atoms do, but hybridisation between Pd 
Ad and N p orbitals leads to formation of ferromagnetic 
order with Pd as the main contributor (not N) to the 
total magnetic moment [2ll |. 

Demonstrated here ferromagnetic order is always more 
energetically stable than antiferromagnetic and paramag- 
netic state, what makes these materials possible candi- 
dates for spin injection in spintronic devices. Calcula- 
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tions of the total energy indicate that this class of ma- 
terials can exist in stable or mctastable phase. Their 
highly interesting magnetic properties should encourage 
experimentalists to stabilize these materials in properly 



coordinated structures via vacuum molecular beam epi- 
taxy, chemical transport or laser deposition on suitable 
substrates. 
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LiN (NaCl, ZB & WZ) E vs a NaN (NaCl, ZB & WZ) E vs a 
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FIG. 1: Energy E (in Ry) vs. lattice constant a (in A) for 
LiN, NaN, KN and RbN in three crystal structures (squares 
- RS, pentagons - WZ, diamonds - ZB, all for ferromagnetic 
phase; open symbols - the same, but for paramagnetic phase) . 
Plots for paramagnetic energies are similar, therefore we give 
them only in the case of LiN. 



ao(A) ^tot(MB) Mi-^(mb) (J.n(hb) Mint(MB) (Ry) 

LiN 



(RS) 


4.36 


2.00 


0.06 


1.72 


0.21 


-0.025 


(WZ) 


3.31 


2.00 


0.03 


1.17 


0.80 


-0.102 


(ZB) 


4.68 


2.00 


0.04 


1.51 


0.44 


-0.052 


NaN 


(RS) 


5.03 


2.00 


0.01 


1.73 


0.26 


-0.061 


(WZ) 


3.91 


2.00 


0.02 


1.47 


0.51 


-0.141 


(ZB) 


5.49 


2.00 


0.01 


1.67 


0.31 


-0.074 


KN 


(RS) 


5.82 


2.00 


0.02 


1.85 


0.12 


-0.075 


(WZ) 


4.51 


2.00 


0.03 


1.75 


0.22 


-0.160 


(ZB) 


6.38 


2.00 


0.03 


1.85 


0.11 


-0.080 


RbN 


(RS) 


6.12 


2.00 


0.03 


1.83 


0.13 


-0.074 


(WZ) 


4.75 


2.00 


0.03 


1.73 


0.24 


-0.156 


(ZB) 


6.70 


2.00 


0.03 


1.83 


0.13 


-0.077 



TABLE I: Equilibrium lattice constants ao (in A) together 
with total /itot , partial atomic-site resolved fijA , /in and inter- 
stitial fii n t magnetic moments (in /xb) as well as total energy 
difference AE t ~^ t (Ry) between ferro- and paramagnetic state 
(in Ry) for RS, WZ and ZB crystal structure. 





K P (eV) 


ET (eV) 


Er (eV) 


LiN (RS) 


8.05 


6.26 


0.2 


LiN (WZ) 


6.56 


4.53 


0.9 


LiN (ZB) 


6.55 


4.52 


1.2 


NaN (RS) 


5.17 


3.12 


1.0 


NaN (WZ) 


4.38 


2.20 


1.7 


NaN (ZB) 


4.44 


2.37 


1.6 


KN (RS) 


4.54 


2.37 


2.0 


KN (WZ) 


3.81 


1.72 


2.0 


KN (ZB) 


3.77 


1.77 


2.0 


RbN (RS) 


4.22 


2.06 


2.0 


RbN (WZ) 


3.48 


1.42 


2.0 


RbN (ZB) 


3.45 


1.39 


2.0 



TABLE II: Energy bandgaps for majority up-spin (-Eg P ) and 
minority down-spin (Eg) states together with half-metallic 
bandgaps (Ef M ) for three crystal structures. 



LiN (NaCl, ZB & WZ) E(MTOT) 



NaN (NaCl, ZB & WZ) E(MTOT) 




1 1.5 
MMTOT [u B ] 



cc 

LJJ 




1 1.5 
MMTOT [n B ] 



KN (NaCl, ZB & WZ) E(MTOT) 



RbN (NaCl, ZB & WZ) E(MTOT) 




cc 

LU 




MMTOT [|i B ] 



MMTOT [n B ] 



FIG. 2: Total energy E (in Ry), relative to the minimum 
energy, for LiN, NaN, KN and RbN in equilibrium for three 
crystal structures (squares - RS, pentagons - WZ, diamonds 
- ZB) as a function of the assumed net magnetic moment fi to t 
(in hb )• Formation of integer magnetic moment p to t = 2/jb 
stabilizes the crystal structure by increasing chemical bonds 
between I A and N atoms. 
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LiN (NaCl, ZB & WZ) MTOT vs a 



NaN (NaCl, ZB & WZ) MTOT vs a 



3. 

h- 
O 




KN (NaCl, ZB & WZ) MTOT vs a 



RbN (NaCl, ZB & WZ) MTOT vs a 




FIG. 3: Magnetic moment /i to t (in /ib) vs. lattice constant 
a (in A) for LiN, NaN, KN and RbN in three crystal struc- 
tures (squares - RS, pentagons - WZ, diamonds - ZB). The 
lines are to guide the eye. Critical values of lattice parameter 
a for compression, below which the total magnetic moment 
/itot = 2/iB disappears, can be seen for each compound in 
three crystal structures. For example, in the case of LiN the 
critical values of a are 2.75 A, 3.75 A and 4.25 A for WZ, ZB 
and RS structures, respectively. 
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RbN RS minority spin 
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RbN WZ majority spin 
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FIG. 4: Spin resolved electron energy bands for RbN in differ- 
ent crystal structures. High symmetry points from the Bril- 
louin zone for two crystal structures are indicated. High- 
est occupied majority spin (up) subbands, shown in the left- 
hand-side diagrams, are separated from the Fermi level Ef by 
half-metallic bandgap E^ M = 2eV. The values of bandgaps 
separating occupied (or partially occupied) and empty states 
for both majority E^ p and minority (down) E^ n spin sub- 
bands (shown in the right-hand side diagrams) are listed in 
Table II. Corresponding total and partial DOS are presented 
in Figure 6. 
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FIG. 5: Spin resolved electron energy bands for RbN and 
NaN in ZB crystal structure. High symmetry points from the 
Brillouin zone for ZB crystal structure are indicated. High- 
est occupied majority spin (up) subbands, shown in the left- 
hand-side diagrams, are separated from the Fermi level Ef by 
half-metallic bandgap E^ M = 2eV. The values of bandgaps 
separating occupied (or partially occupied) and empty states 
for both majority E^ p and minority (down) Eg n spin sub- 
bands (shown in the right-hand side diagrams) are listed in 
Table II. Corresponding total and partial DOS are presented 
in Figure 6. 



11 



DoS - RbN (NaCl) 



DoS - RbN (WZ) 



e/3 
o 
Q 



0.5 
0.25 


-0.25 
-0.5 
5 
2.5 


-2.5 
-5 

0.05 


-0.05 

7.5 
5 



-5 
-7.5 



-5 



-2.5 



■ Rbe„ 

- Rbt 2 ° 










Rbp " 

N p 


1 ^7" 




Rb s 








■ ^ 


total DoS ! 







2.5 5 
E[eV] 



7.5 



10 



12.5 



CO 
o 
D 



0.5 
0.25 


-0.25 
-0.5 
5 
2.5 


-2.5 
-5 

0.05 


-0.05 

7.5 
5 

2.5 


-2.5 
-5 
-7.5 



: a „ 


p — ■ ; v 


Rb eg 

. Rbt2g 






Rbp " 

Np 






— u 




Rbs / 

Ns If 




: 1 




total DoS - 





-5 -2.5 2.5 5 7.5 10 12.5 
E[eV] 



DoS - RbN (ZB) 



DoS - NaN (ZB) 



0.5 
0.25 


-0.25 
-0.5 
7.5 
5 
2.5 


-2.5 
-5 
-7.5 

0.05 



-0.05 

7.5 
5 
2.5 

5 -2% 
-5 
-7.5 



CO 



eo 



-2.5 







- Rbeg 

. Rbt2g 




1 


Rbp I 

N p - 






— L_ 




Rbs \ 

Ns ! 


j — * v- y ^ 




total DoS _ 







2.5 5 
E[eV] 



7.5 



10 12.5 



c/) 
o 
Q 



0.05 


-0.05 

7.5 
5 

2.5 


-2.5 
-5 
-7.5 

0.05 



-0.05 

7.5 
5 

2.5 


-2.5 
-5 
-7.5 





Rbeg 

. Rbj2g ■■ 




—T - .... A ly L.I ^ 






Rbp — : 

Np - 










Rbs y 

" Ns » 








total DoS . 


Y 


^ 



-2.5 2.5 5 7.5 10 12.5 
E[eV] 



FIG. 6: Total and partial site and symmetry projected DOS 
for RbN and NaN in different crystalline structures. Angular 
electron density distribution in the real space for the d x i_ y i 
state is in x and y axis, for the state d iz i_ r i in z axis, what 
increases the overlap with N p electron states in the RS struc- 
ture and explains the hybridization between N p and e 3 states. 
On the other hand, angular dependence of t^g states is fa- 
vored in the ZB structure because of tetrahedral surrounding 
of neighboring atoms and supports hybridization between N p 
and tig states. The hexagonal z axis in WZ structure distin- 
guishes both t2 g and e g Rb states as partners for hybridization 
with N p states. 



